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ABSTRACT 


An experimental investigation was undertaken to evaluate 
the utility of the time series modelling technique proposed by 
Box and Jenkins. In this method input - output time series data and 
their correlation functions are used to develop transfer function 
models of the process and also of disturbances affecting the process. 

Simple dynamic, stochastic models were developed for the 
evaporator at the University of Alberta using experimental and simulation 
data. Experimental data was collected using an IBM-1800 process control 
computer. The effects of feed flow disturbances and the level control 
constants on both models was studies. Identification was carried out using 
both open loop and closed loop data. The results of the identification study 
demonstrated that almost identical transfer function models can be 
obtained from open and closed loop data. 

Using these dynamic stochastic time series models, minimum 
mean square error (MMSE) controllers for product concentration were 
designed. The performance of these controllers was compared with a 
well tuned proportional - integral (PI) controller using both 
stochastic and deterministic disturbances. The simulation study indicated 
that MMSE controllers generally produce smaller sum of squares of 
errors for both types of disturbances. For the experimental study the 
gains of the controllers had to be reduced to produce satisfactory 
control. Good concentration control was produced in spite of severe 
periodic vacuum disturbances. The stochastic controllers based on 


experimental closed loop data produced unstable performance, however 
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the MMSE controllers designed using closed loop simulation data 
performed very well. The simulation study also indicates that MMSE 
control can be successfully applied to the case of setpoint 


changes. 
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CHAPTER ONE 


INTRODUCTION 


The design of a process control system necessarily involves 
understanding the dynamic behaviour of the process. The widely 
used classical control theory ignores the stochastic nature of the 
disturbances involved in chemical processes. A mathematical model 
for such processes that can be related to real data must take into 
account the dynamic relationship between process inputs and outputs 
and also the noise affecting the system. such models can be 
obtained by combining a deterministic transfer function model with 
a stochastic noise model. 

Of the many methods of stochastic model building the 
"Time Series" technique of Box and Jenkins [1] has gained consider- 
able favour. They use single input - single output, discrete, 
linear, transfer function models to represent the process and 
employ a class of statistical time series models known as 
autoregressive - integrated - moving average (ARIMA) models to 
characterise the stochastic disturbances in the system. These models 
are obtained from process input - Sasa data using iterative, 
non-linear least squares parameter fitting. 

Box and Jenkins have shown that a controller which cancels 
the effect of the forecasted disturbances also minimizes the 
variance of output deviations from the setpoint (i.e., a minimum mean 
Square error controller). The time series model is used to obtain 


an estimate of future disturbances. 
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This thesis presents an application of the time series 
modelling technique to model the double effect evaporator in the 
Department of Chemical Engineering at the University of Alberta 
and to design and implement minimum mean square error controllers. 
This application involves extensive digital computer simulations 


and experimental testing on the actual evaporator. 


1.1 Objectives of the Study 


The objectives of this study are twofold. Although there have 
been several applications of the Box and Jenkins [1] technique of 
model building to model various processes [2 - 6], none of these 
include a detailed evaluation of the effectiveness of this technique. 
Thus the first objective of this thesis is to extensively evaluate 
the utility of this model building technique to model interacting, 
multivariable processes such as the pilot plant, double effect 
evaporator. 

The second objective is to evaluate experimentally the 
performance of the minimum mean square controller for a wide variety 


of conditions including various types of disturbances. 


1.2 Structure of the Thesis 

The thesis is divided into three parts: literature survey, 
process identification and controller evaluation. 

A brief survey of the literature on time series modelling 
depicting the theoretical development and applications of this 
technique to various engineering processes is presented in Chapter 


Two. 
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Chapter Three contains the results of model identifica- 
tion and estimation for the evaporator system using both simulated 
and experimental data. The evaporator system was perturbed under 
various open loop and closed loop conditions to study the dependence 
of estimated models on these conditions. 

In Chapter Four the performance of controllers designed 
using these models is evaluated using various types of disturbances. 
Both experimental and simulation results are presented. 


The conclusions from this investigation are presented in 


Chapter Five. 


CHAPTER TWO 


A REVIEW OF TIME SERIES ANALYSIS 


2.1 (Introduction 

Since Box and Jenkins [1] presented their approach to 
stochastic time series analysis in 1970, several researchers have 
investigated this technique for the purpose of identification and 
control of industrial processes [2 - 7]. In this chapter the 
Box-Jenkins technique is briefly described and the important steps 
and equations are presented. Modifications that have widened the 
applicability of this technique are also reetiaed 

Although time series analysis finds extensive applications 
in economics and other fields for forecasting purposes, this survey 
is restricted to applications involving the identification and 


control of chemical processes. 


2.2 The Technique of Time Series Analysis 


In 1970 Box and Jenkins [1] published a new approach to 
linear stochastic modelling using statistical time series maa vote 
The time series considered here are a sequence of values at equally 
spaced points in time. Box and Jenkins have presented a very 
straightforward technique.for obtaining discrete dynamic process 
models from the input/output data. This technique consists of a 
three stage iterative procedure [7]: 

the Identification 
2 Estimation or fitting 
3% Diagnostic checking 


In the identification stage several tentative model forms are 
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obtained from the correlation patterns of the data and prior 
knowledge of the process. The data is then fitted using 
efficient methods assuming one such model form to be true. The 
residual sequence thus obtained is analysed to detect any lack 
of fit as well as its possible cause and remedy. 

Discrete transfer function models are used to describe 
the process dynamics and a class of time series models known as 
autoregressive ~- integrated - moving average models (ARIMA) AS 
used to characterise the noise. 


The input/output relation can be written as, 


t t= Ze t=2 r t-1 o §6t-b 
faa * . ay Se . 
4) ene ar ovenes nh (2.1) 
where, 

x. = value of output deviation from steady state 

value at time t, 
X = value of input deviation from steady state at time t, 
N. = deviations in the output due to noise, 
a. = zero mean white noise, 
b = time delay expressed as an integer multiple of the 


sampling interval. 
oe are autoregressive parameters. 
; fw, } are moving average parameters. 
Equation (2.1) can be written in operator notation as, 
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Suppose that the noise signal Ny is represented at the output 


of a dynamic system whose input is white noise, a,: 


8 (B) 


N 


. a (2.3) 
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Then combining eqns. (2.2) and (2.3) gives the general form of 


the model, 
b 
w (B)B 0 (B) 
. Se Ba ee eA C224) 
t 6 (B) re 4 (zy Vo t 
Pp 
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V = the difference operator (e.g. VX =X - X mg 


V = ), the summation operator. 
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In eqn, (2.4) 6 (B)w, (B)B” represents the process transfer 
function model of orders (r, s, b), whereas p, (BO, (BY 
characterizes the disturbance model, Neo of model orders 

Qa hy Oe Ny is the consolidated effect of noise from all 

sources such as other input variables and disturbances. It also 
represents the effect of modelling errors. Ny and ay are 
assumed to be uncorrelated with Xi Equation (2.4) is represented 
in block diagram form in Fig. 2.1. It is clear from Fig. 2.1 that 
N. represents the value of Y,. when there is no control 

Giechmat t=, 0). 


(process model) 


FIGURE 2.1 OPEN LOOP PROCESS MODEL 
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2.2.1 Model Identification 


From Fig. 2.1, it follows that 


‘ w 6B) BP 
ia Eon (es xX, (225) 


KG 


Now Y. can also be represented in a linear filter form as: 
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= V(B) X,. (2.6) 


The linear operator V(B) is the transfer function and the 
weights ve: vy eee. are called the impulse response function of 


the system. 


Rearranging eqn. (2.5) gives: 
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By comparing eqns. (2.6) and (2.8), the following identity is 


obtained [l, § 10.2.2] 


Veen eee: Sk 
( 6B ) ive + Wy Betaree se) 


We = 0 a. oS Ae 

vy = rT Mees “= 55 ve eae 4 Se ee ae W, j =b 

VE = Sy ee ae 5, NE eee eet oF eae - en j= be, 

Me = 8) ae ete 6, Vey “ie ar ey fee Deets 
C750) 


Equation set (2.10) indicates that if the transfer function 
model is of orders (r, s, b), then the impulse response weights We 
consist of: 


1. The first b weights Vs are zero. 


2. The next s-r+l values do not follow any fixed 
pattern. 


3. The values of Ne Corns betes to 1s etol Low 


the pattern dictated by the a order difference 
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equation which has y starting values, 
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To obtain consistent parameter estimates it is desirable to 
have the input X's uncorrelated with N,'S. This can be achieved 
by having an independent random process generating the input signal 
Xe When Xx is not a white noise sequence computational 
simplicity results by "prewhitening". If an appropriate model for 


xX, is available, then it can be used to transform the auto- 


correlated input series Xi» into a white noise series, a, , 
-] £ 
Oe (B) 8. (B) Xx = a | (CAL) 


where the model for Xx is 
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When such a prewhitened input is used, the ue weights are 


proportional to the cross correlation, p (K) [l, p. 380]. 


aB 


p_,(k)o 

ig Sot Ms a Ne hoo ke hte (2.12) 

J (G) 

(OY 

where 
Bil COne ni) (Be) a) 
t a t 6B 

= Ze 
Pg 6k) 5 ( ) 


E(a, - on E (8 - Up) 


11 


t=1 * 

N 

BAN y, OURRe eS ") 

%% a E(o, Uy) TOES z (a, - My) 
t=1 

and 
E 2 : 
peyn con (B) 0. (B) on (2514) 


Knowing the ae a tentative model form can be identified 
and initial estimates of the (r + s + 1) parameters can be obtained 
From the first (r + s + 1) equations in eqn. (2.10). 

The identification of the noise model is done using the 
auto and partial correlation functions of the residuals, (N. a A ee KORA 
When the residuals indicate nonstationary behaviour (i.e. an 
apparently time varying mean) by differencing the residual series 


q 
one or more times stationary series, Nes can be obtained 


where, d = eeres of differencing. 

Auto and partial correlations for this stationary series 
are then used to obtain a noise model. 

A detailed discussion of these identification procedures 


can be found in Box and Jenkins [1l, pp. 174-194 and 377-388] . 
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2.2.2 Retimation of Parameters 
The estimation of parameters is efficiently carried out 


by minimizing the conditional sum of squares of the residuals. 


S(wo» 55> $5> 9) = min S(w, 6, $, 8) (2.16) 


where 


n 
S(w, 8, OF 6) = acre ior W, o, 6 |x, ves a.) C2 a) 
n = number of points in the series 


tn A Oe starting values of the series re 


(@) ic 


If the first (u + p) values of a, are assumed to be equal to zero, 


where u = max (b +s, rT), then eqn. (2.17a) can be written as 


S(w, 8, 4,9) = = ar(b, sw, 4, 8) (2.17) 
t=utptl 
This maximizes the likelihood function of the parameters. Numerical 
optimization techniques are used to obtain parameter values that 


minimize the residual sum of squares. 


2.2.3 Diagnostic Checking 


The parameter fit is quite good if there is no significant 
cross correlation between the prewhitened input sequence fa, } and 
the residuals fat. If the auto and partial correlations of the 
residuals are within the confidence limits, then the noise model 


is acceptable. Two quantities ae and Qe which are defined 
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below are approximately distributed as y7 with f£ degrees 
of freedom. They are used to test model adequacy by means of a 


chi-square test. 


Abe 
Svalf) = ine. Pua (®) (2718) 
k=o 
é 2 
Q,(£) mT eed Yee tk) C2519) 
k=] 
where, 
n = total number of data points. 
K = number of lags for which the values are calculated. 
f= (K +41) - (r+S +41) for Se 
and foee(h) vero + ade for Q, 
Sa and Q. are compared with values in a table of 


percentage points of ve to obtain an approximate test of the 


hypothesis of model adequacy. A discussion of these tests can be 


Loundeane Box and Jenkins i> lop. »392,-,.395, 522°): 


2.2.4 Closed-Loop Process Model Identification 


Often process operating data is used in the model 
identification study. In many cases it is difficult to obtain 
open loop operating data because open loop conditions may lead to 
undesirably large deviations and an inferior quality product. 


Also from a stability point of view, it may not be desirable to 
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operate the plant under open loop conditions. 

The problem of identification under closed-loop conditions 
has been considered by Box and MacGregor [8, 9] and Graupe [10]. 
They have shown that to identify an unknown system it is often 
desirable to add an artificially generated signal De to the input, 


hoe einem dither tsiena lly: is a stochastic sequence with zero 


iE tC 


mean and is uncorrelated with the error, eny 


From fig. 2.2, the feedback controller equation is, 
X = C(B)e, + D. (27420) 
where C(B) is the controller transfer function and 
Orta Cr V(B) X, + y(B)a, (2023) 


t 


Tnweqnan(272)) v7 is the deviation of output from 


steady state. 


process 


V(B) 


controller 


C(B) 


FIGURE 2.2 CLOSED LOOP SYSTEM WITH DITHER SIGNAL, D 
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Using eqns. (2.20) and (2.21) and substituting for xX, gives 
1 - V(B)C(B) e. = (B) ay + V(B) D,. Cee) 


For the case where D, = 0 and the model orders are 
not known, then unambiguous identification will not in general be 
possible. If either V(B) or w(B) is known, then using auto- 
correlations of e, it is possible to arrive at a correct model 
for the other one. If the correlation among parameters is high, 
difficulties can arise which can be greatly simplified by using 
a non-zero De. signal [8]. 

When D, + 0, eqn. (2.22) can be rearranged as 


0 UO y(B) 
eT he MEDENTCOOED ae ee 


Comparing Fig. 2.3 with Fig. 2.1 it can be seen that for 


non-zero JD a procedure similar to open loop process model 


t’ 
identification can be used where the dither signal De. is treated 
as an input. Here instead of relying on autocorrelation 

properties of ere the cross correlations between dD. and om 


and the autocorrelation of e, can be used for parameter estimation. 


16 


VAG) 
1 = V(B) C(B) 


V(B) 
1 = V(B) C(B) 


FIGURE 2.3 REPRESENTATION OF CLOSED LOOP PROCESS IN OPEN LOOP FORM. 


Hong and MacGregor [5] recommend the use of D-optimal 
dither signals (signals that minimise the volume of joint confidence 
region of the parameters) for obtaining better parameter estimates. 
When the variance of the input or output is constrained, first order 
autoregressive dither signals with appropriate values of the 


parameter, 06, are shown to be D-optimal. 


2-3 Design of a Feedback Controller 


Once the noise and transfer function models are obtained, 
the next step is the design of a controller to maintain the output 
in proximity of the setpoint value. Box and Jenkins [l, 8 12.2] 
have suggested a procedure to design a controller that minimises the 


mean square error between output and setpoint. 
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From Fig. 2.1 it can be seen that for constant input 
conditions the output deviations in the error are due to disturbance 
term Nee Therefore, Nye represents the total effect of 
disturbances and the total effect of compensation is the term 


5! (By w(B) xee 


If it were possible to set, 


es 
Xia = <(B) w § (B) Nap ) G2e24) 


where sar denotes the calculated control action after the error measure- 


ment, e, is obtained, the effect of disturbance would be cancelled, but 


since b is positive, the value of Neon isnot available at time t. 


However by using an estimated or forecast value, Nb), the minimum mean 


square error can be achieved. 
Let, 


(b) = N. - Ny, ()) C225) 
and 


Shige 
K4 = -5(B) w (BY) N,_, (b) (2.26) 


where 


“A 


Ny p>) is a b step ahead forecast of N, at 


time t-b and e, p 6b) a is the measured error. 
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Since 


sia teh 
N. = Pan CB) = -0.C 8) ay 627233) 


N. can be written in terms of past, present and future values 


of a as, 


Ne = )\(B) ab + n(B) ay. G27) 
Then, 
N,.(b) = n(B) ay 
e, (db) = 4(B) a, (2.28) 
From eqn. (2.28) 
‘ Pan(B) ve _ n(B) | 
N, (b) “55(3) ep 6b) Sein) e, (2.29) 


The previous derivation is subject to certain restrictions. 


It is possible to set 


a = -§(B) w (B) Neb (2.24) 
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mee om. (B) 
w(B) d(B) 


FIGURE 2.4 MINIMUM MEAN SQUARE ERROR (MMSE) FEEDBACK CONTROL SYSTEM 


only when the roots of 6(B) lie outside the unit circle, or, in 
other words, when 46(B) is invertible. Similarly the noise model 
needs to be invertible, with the roots of ¢(B) lying outside the 
unit circle. 

Further analysis of the minimum mean square error controller 
indicates that in some cases the MMSE controller will produce large 
adjustments in the input Xs This is often the case when the 
polynomial 6(B) has a root near the unit circle. Bacon and Howe 


[3] note that in some cases the MMSE controller can produce unstable 
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system performance even when the process model is known. If 

large variations in the input cannot be tolerated, redesigning the 
controller is necessary such that the input variations are 
constrained at the cost of small increases in the output variance 


lesiplls uel as 


2.4 Applications 
The Box and Jenkins [1] technique is essentially a "black box" 


approach; that is, knowledge of the physical process is not necessary 
for model building. Normal plant operating data can be used to 
obtain dynamic and noise models. Because of this simplicity of 
analysis, this technique has been found quite attractive. Several 
investigations have been carried out to verify the applicability of 
this procedure to modelling various processes. 

Tee and Wu [2] applied this technique to a paper making 
process. The input variable was the scale reading an the stock 
gate opening and the output was the basis weight erin paper 
produced. By keeping a constant gate opening, they developed a 
disturbance model based on measurements of the basis weight at 
20 minute intervals. Later they obtained a combined dynamic and 
disturbance model, using a experienced operator to perturb the 
input variable and assuming open loop conditions. A MMSE controller 
was designed based on this model and it provided improved control 
over the open loop case. 

Restrepo et. al. [13] used this technique to test the 
feasibility of using time series analysis to model a high order, 


non-linear system over a wide range of operating conditions. 
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The system considered was a dynamic model consisting of 20 linear 
and 10 nonlinear differential equations for a pilot plant scale 
distillation column with eight sieve trays. This theoretical 
model was based on physical knowledge of the system. Input-output 
data was generated from this model to develop transfer function 
models using the Box and Jenkins technique. The step and impulse 
responses of these models were found to be extremely close to those 
obtained using the higher order physical model. 

Bacon and Howe [3] assessed the utility of the Box and 
Jenkins time series modelling technique to model an industrial 
falling film evaporator system and used the resulting model 
to evaluate and modify the existing control system. They found this 
technique adequate to model a full scale chemical process unit in 
terms of dynamic as well as stochastic relationships. The models 
thus obtained were very useful in deriving novel control functions 
that supplement the existing control schemes. However the new 
control technique was not applied. 

The time series technique was used by Wright 
and Bacon [4] to identify the dynamic behaviour of a laboratory 
scale heat exchanger network. This system was excited by perturbing 
the steam control valve position at various inlet water flow rates. 
The effects of changes in process conditions and process nonlinear- 
ities on the estimated models and parameters were studied. Wright 
and Bacon suggest that experience and basic knowledge of the process 
must be used to eliminate models which although satisfactory from a 


Statistical point of view are not adequate for physical reasons. 
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Two applications of model building from closed-loop 
operating data have been cited by Box and MacGregor [8]. The first one 
is a polymer viscosity control system. The other an example 
of a paper making process [2]. Here in an attempt to perturb the 
system in a random fashion, the experienced operator regulated the 
stock gate opening, i.e. he acted as a feedback controller [8, 14, 15]. 
Hence the operating data was analysed once more to obtain models using 
closed loop methods. 

Hong and MacGregor [5] have used the closed loop analysis 
to model a continuous stirred tank heating process operating under 
a discrete PI controller. They note that in order to develop 
confidence in the identified models, the choice of input signal is 
important. | 

The sampling interval [16] can also be chosen appropriately 
to obtain maximum information about the dynamic stochastic behaviour 
of the process. 

Hong and MacGregor [5] and MacGregor et. al [7] have 
designed stochastic and PI controllers using the identified 
stochastic models. In an experimental evaluation using stochastic 
disturbances, these controllers performed better than PI controllers 
designed using conventional tuning methods. 

MacGregor and Tidwell [17] have considered the problem of 
stochastic control where there are constraints on input variations. 
They have discussed methods of designing discrete stochastic controllers, 
+ } xa LOG 


: Berd 
which minimise the quadratic cost function, E Ye 


linear processes under stochastic disturbances. Using these methods 
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they designed a controller for the control of viscosity in an 
industrial polymerization reactor. This controller was found to 
perform satisfactorily. 


A discussion of minimum variance controllers can also 


be found in [18]. 
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CHAPTER THREE 


PROCESS MODEL IDENTIFICATION 


Seen roduct ion 

This chapter presents the various time series models 
obtained for the double effect evaporator using both experimental 
and simulation data. A brief description of the evaporator system 
and the computer programs used is also presented. 

The model identification is discussed in two main parts. 
The first part (§ 3.4 - 3.6) deals with identification under open 
loop operating conditions and the second part (§ 3.7 -3.8) deals 
with identification under closed loop operating conditions. Both 
parts include simulation and experimental results. Finally, 
comparisons of process models using open and closed loop operating 


data are presented in sections 3.7 and 3.8. 


3.2 Description of the Evaporator System 


Figure 3.1 shows a schematic diagram of the double effect 
evaporator in the forward feed mode of operation and a conventional 
multi-loop control system. The eight inch diameter, calandria type, 
first effect contains thirty-two, QO.75in.0O.D. tubes and operates with 
a nominal feed rate of 4.5 lb/min of 3 percent aqueous triethylene 
glycol. It is heated by 1.8 lb/min of saturated steam. The first 
effect overhead vapour is utilized to heat the first effect bottoms 
which flows under forced circulation in the three, 6 feet long, one 
inch 0.D. tubes of the second effect. The second effect operates 


under vacuum. The product is about 10 percent aqueous triethylene 
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glycol. 

An IBM 1800 real time digital computer is interfaced to 
the evaporator. Conventional multi~loop control of the evaporator 
is achieved using the IBM 1800 Direct Digital Control (DDC) package 
and single and cascaded control loops. Although product concentration, 
C2, is the most important controlled variable, the liquid levels, 
Wl and W2, must also be controlled within operating levels. For 
the present study the liquid levels were controlled by manipulating 
bottom flow rates, Bl and B2, respectively. The product 
concentration C2, was controlled by manipulating the flow rate of 
steam S, into the first effect. 

A detailed description of the evaporator system can be 


found elsewhere [19]. 


3.2.1 State Space Model of the Evaporator 


A number of process models ranging from tenth order 
nonlinear state space models to the first order transfer function 
models have been derived for the evaporator by previous investigators. 
For the present simulation study the discrete, linear fifth order 
state space model developed by Wilson [20] was used. 


This model can be represented as 


(Let) moe (tf jeereni(t) -t.0d(t) (oe) 


y(t) = C x(t) (San) 


where 


a (59x 1)'state vector 


ix 
ee 
6) 


u is a (3 X 1) control vector 


|a 
}- 
n 


a (3 X 1) disturbance vector 
Voie ee Qjecul) output vector 
t is the sampling instant 

and o, A 


8, C, are constant coefficient matrices 


of appropriate dimensions. 


The variables in eqns. (3.1) and (3.2) are expressed in 
normalized perturbation form: e.g. 


Wl - Wl 
en ese 


1 Wl (Bins) 
SS 


where subscript ss denotes the nominal steady state value. The 


variables in eqns. (3.1) and (3.2) are defined as follows: 


x! = [wi', cl’, Hl, W2', C2" 
uw = S4, eB lise 20 
T 


ge ecb aCh ieee!) 


and 


a =[WI1', W2" C2") 


The nominal steady state values and coefficient matrices 


are presented in Appendix A. 
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The proportional control law used in the closed loop 


tdentification study can be represented as, 


u(t) = K x(t) (3.5) 
where 


K its a 3 X 5 matrix of controller constants. 


3.3 Computer Program for the Simulation Study 


The discrete fifth order state space model [20] was 
simulated on the University's Amdah1-470V/6 digital computer. A 
program, called SIMUL, was written to simulate the response of the 
evaporator for various types of disturbances under conventional 
multi-loop and minimum mean square error (MMSE) controllers. 
Stochastic disturbances such as Gaussian white noise and first order, 
autoregressive signals and step changes in load variables were 
used to evaluate the controllers, 

All the relevant state,control and disturbance variables 
were printed out and punched on cards Pence in plotting programs, 


RBNO1 and RBNO2 [21]. 


3.4 Open Loop Identification Studies 


For the paretor system, the most important output 
variable, C2, is controlled by adjusting the steam flow rate to 
the first effect. Hence process model identification studies were 
carried out between input variable, S, and output variable, (C2. 
Both process and noise models were developed for a variety of process 


conditions. Experimental and simulation runs were carried out for 


two cases: (i) steam perturbations only, and (ii) simultaneous 
perturbations of steam flow and feed flow. Thus the feed flow 
perturbations introduced "noise" into the system. 

The perturbation signals consisted of either Gaussian 
white noise or first order autoregressive noise sequences with 
zero means and various standard deviations. 

The initial model orders and the initial parameter 
estimates were obtained from the cross correlations between the input 
and the output and also from the auto and partial correlation 
functions of the prewhitened series and the residuals. These 
initial parameter estimates and model orders were then used in a 
nonlinear, iterative least squares program, crt, that minimizes 
the sum of squares of the residuals. The stopping criterion employed 
in the estimation program was that the change in the sum of squares 
of the residuals was less than 10mee or the change in each parameter 
value of less than 10h. whichever happened first. 95 percent 


confidence limits for the parameters were also calculated. 


3.5 Identification of Evaporator Models Using Open Loop Simulation 


Data 


Time series models were obtained from data generated from 


the fifth order state space evaporator model presented in Appendix A. 


Each run consisted of 320 data points obtained at 64 second intervals. 


A program made available by Dr. J.D. Wright, McMaster University, 
Hamilton, Ontario. 
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Since these points covered a time period of six hours, the effect 
of slower evaporator modes (i.e. large time constants) were fully 
realized in the input-output series. A large number of data points 
was used in order to produce small confidence intervals for the 
parameter estimates. Also the confidence limits on the correlations 
are narrower (the confidence limits are + 2 + (number of data 
podmee)@ Ne 

For purpose of comparison, a theoretical hee transfer 
function model was calculated from the discrete fifth order model 
using the GEMSCOPE program [22]. To check how closely the estimn- 
ated models agreed with this theoretical model, in some runs para~ 
meter estimation was performed using transfer function model orders 
which were the same as those for the theoretical model. The effect 
of different model structures and different input signals on the 
parameter estimates was also studies. 

Table 3.1 lists the run conditions for open loop 
simulation runs SS-1 to SS-6, including standard deviations of the 


input sequence, the first order autoregressive filter 


Oo 

input’ 
constant, dg» and the actual input and output standard deviations. 
Figure 3.2 presents a base case simulation run, #SS-4, where the 


two levels are under proportional control using the control law, 


uk) E= KX x(k) , where 


u, | (3.6) 


0 0 


le 
U1 
ict 
oo 
pas 
oo 
is) 
os 
ae 
Ul 
em | 
e 


and 


30 


ad 
ei) i 
Is Fe ey 4 ‘ 
: a 
Mt 1 f 
’ 


eae? 


i and seadiinidl t 4% a: 


TABLE 3.1 


CONDITIONS FOR OPEN LOOP SIMULATION RUNS 


(Steam disturbances, evel control only) 


Steam (1b/min) 


Concentration 


oH 


32° 


Vu AL 


Ww 


Ce? 


430 Ni 


OQtee 


7-SS Nha dod VLVG NOILVINWIS dOOl NddO Z°€ TNOTA 


SSUNIA NI SAIL 


CCT OE O72 OST 08 ©) 


1 


are 


' 
Ste a iene peri 


— 


—— 


hanphfh he 


| 


nae | 


i: 


| 
| 


4) 


NIWGT NT MOU WVALS 


a 


cr ey Z 


v 


NIWA NI ALE 


age) 


S3IPNIN NI SIL 


or ores _ Ove OST CS G 


OU 6 


TL 


a 


TAT) 


GOH «NI NI 


DOAN Jd NI NONT) 


In this series of runs, only steam flow rate was 
perturbed. Table 3.2 compares the parameter estimates obtained using 
dimensionless x and y for these runs with the theoretical model. Table 
3.3 lists parameter estimates and the 95 percent confidence limits for 
each parameter. Appendix D lists the 95 percent confidence limits for 
each parameter and also additional models that were fit with their 
parameter estimates. Only the best fitting models are presented in 
Tables 3.2 and 3.3. 

Note that runs SS-4 to SS-6 have approximately equal 
input standard deviations but different input series structures. 
In Figs. 3.3 and 3.4 typical response data for runs ss-5 and SS-6 
are presented. 

Table 3.4 presents the various run conditions for runs 
SSF-1 to SSF-4 where both steam and feed flow disturbances were 
employed. Typical response data are presented in Figs. 3.5 and 
3.6. The steam and feed flow disturbances were either gaussian or 
autoregressive noise sequences. For these runs, two gaussian random 
number sequences were generated, one for steam and the other for 
feed flow. The actual disturbance signals were then calculated by 
passing these sequences through first order autoregressive filters. 

Run lengths of 320 data points were selected in order to 
produce small confidence intervals for the parameter estimates. 
Since these points cover a time period of six hours (sampling interval 
= 64 secs.), ee effects of the slower evaporator modes (i.e. large 


time constants) are fully realized in the input-output series. 
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ano edeeeChoice of Input Signal 


For model identification it is desirable to have an 
input signal which excites all modes of the process go that the 
information ene ene of the input and output series is large and 
a satisfactory model can be obtained from the data. It is also 
desirable to have input signals that minimise the volume of the 
joint confidence interval of the parameter estimates (i.e. D-optimal 
signals). This objective may require very large input variations 
which are not always possible in practice. For cases where the 
input variance is constrained, Hong and McGregor [5] have suggested 
the use of first order autoregressive input sequences with positive 
parameters(i.e. ¢ > 0). 

In three runs, SS-4 to SS-6, the standard deviation 
of steam, Og, is fairly constant so the effect of the autoregressive 


parameter > on the output standard deviation o can be evaluated. 


CZs 


Table 3.1 indicates that o increases sharply as the autoregressive 


C2 
parameter $ increases and o, remains almost constant. For a value 


OfsnoCe@nO0. 5S. Cnen © value is almost twice as large as the value 


674 
EOTsthnersoe=a008(Orawhite noise) case. “Similarly, for: ¢ = 0.75 

the standard deviation is three times that of the run with $6 = 0. 
From Table 3.3 it is clear that for autoregressive input signals, the 


confidence limits on the parameter estimates become smaller as $6 


increases indicating more confidence in the parameter estimates. 


eres Comparison with Theoretical Process Transfer Function Model 


The parameter estimates and the theoretical model 


parameters are presented in Table 3.2. The 95 percent confidence 
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intervals for the parameters appear in the Appendix D. It can be 
seen that when the assumed model orders are correct (runs SS-1 to 
SS-3), the autoregressive parameters {6,} are exactly the same 


j 


as those of the theoretical model but some of the moving average 


parameters wy} differ significantly ( v 20%) from the theoretical 
values. The time series model produces residuals with a small sum 

of squares of residuals but the cross correlation between prewhitened 
inputs and the residuals are significant, as indicated by the large 
values of oa However, the cross correlation values do not show 
any particular pattern. The S va values for runs SS-1 to SS-3 

are much larger than the value of 25.0 which corresponds to the 

95 percent confidence level for 15 degrees of freedom. This indicates 
little confidence in the models, whereas the sum of squares of 
residuals indicates a good fit in the least square sense. 

In Table 3.2, for runs SS-l and SS-2 the Q. values 
indicate adequate noise models but since the Sa values indicate 
little confidence in the process models and the noise is dependent 
on Piencorredtions of the process models, noise model adequacy is 
of little significance. 

Runs SS-4 to SS-6 with "incorrect" model orders 
indicate more confidence in the models as can be judged from the 


small Sua values; however they produce higher sum of squares 


of residuals. 
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3.5.3 Models from Open Loop Simulation Data 


In the simulation runs of Tables 3.1 and 3.2 where only 
the steam is disturbed, the residuals represented modelling errors. 
In order to have a more meaningful noise term, the feed flow rate 
was also perturbed in a random fashion which produced deviations 
in the output variable that were not related to the input (steam) 
signal. The extent of change in model orders and parameter estimates 
in presence of feed flow perturbations was studied by comparing runs 
SSF-1 to SSF-4 in Table 3.4 with runs SS-4 to SS-6 where feed 
disturbances were absent. 

Figure 3.7a presents a typical plot of cross correlations 
between prewhitened input and transformed output for run SSF-2. Up 
to lag zero there are no significant cross correlations indicating 
a dead time of one sampling interval. From lag k= 1 onwards, the 
cross correlations increase gradually and are almost constant between 
Kee eou COmeKe@= 015, Lt was difficult. to analyse this plot to: get 
initial model orders. Consequently knowledge of the theoretical 
evaporator model indicated that model orders of (r, s, b) = (2, 1, 1) 
would be reasonable. The initial parameter estimates thus obtained 
were used to calculate initial estimates of Nye The values of auto 
and partial auto correlations ene N, indicated a mixed first 
order autoregressive-moving average model. All of the initial 
parameter estimates were used in the nonlinear, recursive, least 
squares, estimation program to calculate the final parameter 
estimates. Figure 3.7b, c and d show the cross correlation, auto and 
partial auto correlations of the residual series. No systematic 


patterns are apparent but there are one or two correlation values that 
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TABLE 3.4 


- 


CONDITIONS FOR OPEN LOOP SIMULATION RUNS 


( Steam and feed flow disturbances, level control only) 


Feed Concentration 


Steam 


(1b/min) (1b/min) (% glycol) 
ES input °F ay ° C2 
0.50 0 0.485 0.65 
OF 5 0neer Om 0'sn5 31) hace 
0.50 O75 Uo S ee 
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FIGURE 3.6 OPEN LOOP SIMULATION DATA FOR RUN SSF-2 
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a. Cross correlation function, Pag 6K) 


c. Autocorrelation function d. Partial autocorrelation 


of residuals,y, (k) function of residuals ,¢, (k) 


FIGURE 3.7 CORRELATIONS FOR SIMULATION RUN SSF-2 
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lie outside the 95 percent confidence limits. In Table 3.5 the 
value of Syl) = 28.87 for run SSIF-2. indicates that the 
model is within the 96 percent confidence limits and similarly 
the value of Q, (18) = 30.08 also indicates 96 percent confidence 
limits, 

For runs both with and without feed flow disturbances 
(i.e. runs SS-4 to SS-6 and SSF-1 to SSF-4) the process 
transfer function models with model orders (r, s, b) = (2, 1, 1) 
were found to be adequate. The 5 ae values for these models were 
generally within 95 percent confidence intervals. With an increase 
in the number of parameters, the fit is generally better (lower chi- 
square values) but at the cost of wider confidence intervals for 
the parameter estimates. 

Pnvsicallyewas. Clem lsat), model has a dead time of one 
sampling interval and means that the current product concentration . 
depends on two past values of steam flow and two past values of 


product concentration: 


tae 
ii 
KS 
+ 
On 
ta 
+ 
S 
ms 
I 
= 
— 
mS 


t ier tel Ze t=e of yet t-2 
gh eek Tet imd-b26 tb 22 (3.8) 
where 

Y. = normalized deviation of product concentration from 
steady state value at time t; 

x, = normalized deviation of steam flow rate from steady 
state value at time t;3 

51> 55> Wo» Wy are constants; 

Ni [noise at tame.) \t. 
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Table 3.6 presents the parameter estimates for run 
SSF-3 assuming different model orders. For different noise models 
the transfer function model parameters Ao are also different. For 


“a 


example, for the second order autoregressive noise models 54 


changes from 1.509 to 1.44 and 54 changes from -0.51 to -0.47. 
On the other hand the ,} are fairly close (within 2 percent of 
each other) except for the runs where no noise models were fitted. 
The variations in parameters indicate their sensitivity to the 
assumed noise model orders. 

The (2, 1, 0) noise model is considered to be the best 
model for this run because it provides small values of chi-square 
statistics (S and Q), (both S and Q are within 95 percent 
confidence intervals) and has one less parameter than the (2, l, 1) 
noise model which produces slightly smaller chi-square values for the 
auto and cross correlations. Other noise models are inadequate 
because they produce high values of Q. For other runs the parameter 
estimates using different noise models and also their confidence 
intervals are presented in Appendix D. 

Runs SSF-3(F) to SSF-3(H) illustrate the effect of 
assuming different process models. For run SSF-3(F), the order of 
S is increased from 1 to 2 and the 6,3 are close to those for 
run SSF-3(A) but the fo.) are completely different. Similarly, 
for runs SSF-3(G) and SSF-3(H) the 65) are close although they 
are quite different from all other runs. For these runs, two of the 


three roots of the polynomial, 6(B), are approximately equal to 


those of run SSF-3(A) but this is not the case for {w(B)}. It 
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can be seen that the process can be adequately described with a 
lower order model as the lower order model also satisfies the 
chi-square criteria. 

By comparing runs SS-4 to SS-6 (Table 3.2) and runs 
SSF-1 to SSF-3 (Table 3.5) it can be observed that the auto- 
regressive model parameter estimates, 64), decrease by up to 
20 percent when feed disturbances are introduced, while the moving 
average parameter estimates, a}, are fairly close. In run SSF-4 
where the feed disturbance magnitude is twice as large as the steam 
flow disturbance, the 64) decrease even furthers Thus it seems 
that as Op/ 5, increases, the 64) decrease but the fa} remain 
fairly constant. 

Looking at the noise models, one sees that the N. series 
is stationary when feed disturbances are absent and the fitted noise 
models have order of the differencing (d) equal to zero. But in 
the presence of feed flow disturbances they exhibit nonstationary 
behaviour and differencing the residual series N, once produces 
a stationary series. As expected, the noise models and parameters 
are completely different for these two types of disturbances. 

For runs SSF-1 and SSF-2 the noise models have the same 
orders and the autoregressive parameters are close. The values of 


n 


moving average parameter 81 have the same signs and are comparable 
in magnitude. For runs SSF-3 and SSF-4 where the input series 


structures are similar ($. and bp constant), the noise models 


have the same structure and the parameter estimates closely agree. 
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3.6 Identification using Open Loop Experimental Data 


As in the simulation study two sets of experimental runs 
were carried out: the first set with steam disturbances only, and 
the second set with additional feed flow disturbances. Tables 3.7 
and 3.8 summarize the run conditions and typical runs are presented 
Dee ees Oo LOmo de. 

In each set of runs, proportional level control was 
employed. Runs S-l to S-4 and SF-1 to SF-4, used controller 
matrix K in eqn. (3.7) which was also used in the simulation 


study. For run SF-5 a different controller matrix, K,, was used: 


Loe O see Oe © 
OFF Os? Oe. 05> £0 


(3.9) 


In the experimental runs 400 data points were collected at 
64 second intervals. For the model identification 320 points 
beginning at the 29th data point were used. The first 28 points were 
not considered in order to eliminate the effect of the initial 
transient caused when disturbances suddenly began. 

The steam perturbations were introduced by changing the 
set point of the steam flow control loop every sampling instant. 
Similarly the total feed. flow rate was changed by changing the set 
points of the feed water and feed solution control loops, so as to 


maintain constant feed concentration. 


Comparison of C2 responses using the simulations (Figs. 3.2 to 
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3.6) and experiments (Figs. 3.8 to 3.12)indicates that unlike the simulated 


responses, experimental responses exhibit time variant means or decaying or 


growing oscillations, indicating different types of process noise. 
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TABLE 3.7 


CONDITIONS FOR OPEN LOOP EXPERIMENTS 


( Steam flow disturbances, level control! only) 


Run No Steam (1b/min) Concentration 


(ZEelycol) 


x 
C2(0) 


* C2(0) is the initial steady state glycol concentration. 


1 Control matrix Ki in eqn. 3.7 was used. 
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FIGURE 3.11 OPEN LOOP EXPERIMENTAL DATA FOR RUN SF-4 
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3.6.1 Process Models from Open Loop Experiments 


Figure 3.13a presents the plot of cross correlations 
between input and output for experimental run S-3. Up to lag zero 
the cross correlation values are well within the 95 percent 
confidence intervals and are almost equal in magnitude. At lag 
k = 1 the cross correlation value increases from 0.047 to 0.077. 
This suggests a transfer function model with b = l. 

As in the case of the simulation data in Fig. 3.7a, the 
cross correlation values gradually increase from lag k= 1 to 
k = 6 and are fairly constant thereafter. Previous knowledge of 
the theoretical evaporator model was used to determine initial model 
orders of (r, s, b) = (2, 1, 1). The estimated noise series Nye 
was calculated using the initial parameter estimates. The auto and 
partial autocorrelations of N, indicated a second order auto- 
regressive structure. After some trial and error, a better fitting 
ARMA noise model with orders (p, d, q) = (2, 0, 1) was obtained. 
The cross, auto and partial autocorrelations of the residual series, 


a are plotted in Fig. 3.13b, c and d. All the correlation values 


t? 
are well within the 95 percent confidence limits and no systematic 
patterns are evident. The one value of 12.99 (from Table 3.9) 
indicates transfer function model adequacy. Similarly, 

Q, 07) = 15.85 indicates noise model adequacy using the chi-square 
test. 

Table 3.9 presents the effect of assuming different model 


orders (process and noise) on the parameter estimates for a typical 


open loop experimental run, S-1l. Runs A-D illustrate the effect of 
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FIGURE 3.13 CORRELATIONS FOR EXPERIMENTAL RUN S-3 
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assuming different noise models. Here a change in the noise model 


has almost no effect on the autoregressive parameters estimates 


{6,}. On the other hand, variations in w, are up to 10 percent 


j 


and Wy 


el y260k 10% 


a 


exhibits very large fluctuations in value from 

3 COnO 02 5X re, For different noise models, the noise 
model parameter estimates have widely varying values but the sum 

of squares of residuals undergoes little change. For runs A-D 

all oe. values are smaller than chi-square value of 27.6, thus 
showing transfer function model adequacy. All the Os values 
indicate confidence limits wider than 95 percent. 

Runs E-H show the effect of using different transfer 
function models. Here the runs have autoregressive parameter. 
estimates that are close to the estimates in runs A-D; however 
the {w,} are quite different. These runs indicate that with more 
parameters (higher orders), the sum of squares of residuals is 
smaller and also smaller Sa are obtained. The use of a noise model 
greatly reduces the sum of squares of residuals. 

Fot open loop runs S-l. to $S-4 .and SF-l to. SF-5, 
the models which provided the best fit are presented in Tables 3.10 and 3.11. 
Additional models and the confidence limits for the parameters are Pr- 
esented in Appendix D. It can be seen that for runs using the level 


COntrod matrix Kh the Ly have remained unchanged with the 


ae 
exception of run SF-2. For run SF-5 using tighter control on 
levels (i.e. control matrix K, in eqn. (3.9)),: the 6) are 


smaller with one of the roots of 6(B) moving away from the unit 


circle, but having little effect on the other root. 
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In the numerator polynomial, w(B), Wo has various 


: COL Oona 1077 for both sets except 


for run SF-2 where it is -0.162 X 10. The value of WF has 


values between 0.61 X 10 


increased by an order of magnitude for the case of simultaneous 


steam and feed disturbances. With the tighter level control in 


run SF-5 the value of OF further increases. Within one set 


of runs the value of Wy exhibits large variations but the order 


of magnitude remains unchanged. This indicates sensitivity of 


n~ 


acre to changes in process conditions. However the autoregressive 


polynomial, 6(B) is sensitive only to the level control matrix. 


The residual noise sequence, N was generally 


t? 
nonstationary, requiring one degree of differencing to get a station- 
ary series. The noise models for all but one of the runs were 
dominated by first order autoregressive structures after differencing. 
once. For run S-3 the noise was best described by a second order 
autoregressive and first order moving average model (model orders 
(2, 0, 1)). In runs SF-2 and SF-4 after differencing once a 
mixed autoregressive moving average holed was observed. 

For the runs with noise model (1, l, 0), the autoregressive 


parameter 1 varied between -0.0968 to -0.226. Other models have 


“A 


a positive moving average parameter 6 with widely varying 


1? 
magnitudes. From this it seem indicated that noise models and their 
parameters are sensitive to changes in process conditions. 

The eniteris used in selecting the best fitting noise 
model was as follows: When two noise models produced nearly equal 


S and Q values, the model having nonzero order of differencing 


(d a 0) was chosen, since the control action of the minimum mean square 
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error controller based on such noise model incorporates integral 
action. The next criterion was to choose a model having the least 
number of parameters from two or more equally adequate models. 

Comparison of models obtained using experimental and 
simulation data shows that the same transfer function model orders 
adequately fit both types of data. The 64) obtain experimentally 
are 10 - 15 percent higher than for simulations, whereas the a) 
are higher in the case of models obtained from simulation study. 


“a lal 


The increase in Wo and w,- are not proportional. The auto- 


1 
regressive parameter estimates, LoVe in the case of the simulation 
change with change in magnitude of feed disturbances, whereas 
they seem to be independent of feed disturbance magnitude in the 
case of models based on experimental data. 

No significant difference is observed between noise model 
orders fitted for experimental and simulation data. Similar Te and 


a8 were obtained in the experimental and simulation data except in cases 


of simulation with steam disturbances only, they were higher. 


3.6.2 Effect of Run Length on the Parameter Estimates 


In order to study the effect of run length on the parameter 
estimates, various esters of data points from run SF-4 were used. 
The results of parameter estimation using process model order (2, 1, 1) 
(0, 0, 0) are presented in Table 3.12. The first four data sets 
contain 120 points each starting at different sampling instants. Sets 
#5 - 8 were similarly chosen with 200 data points in each set. Para- 


meter estimates from these sets are compared with those from the 
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entire run of 320 points (Set #9). 

The parameter estimates in Table 3.12 indicate that the 
| 6) for the runs with 120 or 200 points are approximately 2 - 3 
percent higher than the parameter estimates for the complete run. 
However, the estimates of Wo for the shorter runs were lower than 
the estimates for total run by 50 - 75 percent. The values of a, 
do not exhibit any particular trend, but deviate from the values 
for the complete run. 

The differences in parameter estimates indicate that when 
short run lengths are used, the effects of the slower modes of the 
process are not fully realized and smaller values for the slower 
time constants are obtained. 

Table 3.12 also presents the confidence intervals of the 
correlations for each set. A value of a cross or autocorrelation 
within the interval may be assumed to be equal to zero. If many 
correlation values lie outside these limits, this indicates model 
inadequacy. Runs with larger numbers of data points produce smaller 
confidence intervals (since the confidence interval = 2/(number of 
data points) indicating more confidence in the correlation values. 

In summary, it can be said that runs with large numbers 
of data points are more likely to take into account slower process 
modes and have narrower confidence intervals indicate more 
confidence in the estimated model. Runs with 240 to 400 data points 


would normally be quite adequate for the double effect evaporator. 
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3.7 Identification from Closed Loop Simulation Data 


Process models were also identified from closed loop 
simulation data to evaluate the applicability of time series 
modelling for processes which must be maintained under closed loop 
control. 

As per § 2.2.4 in order to obtain better estimates of the 
model parameters, a dither signal was added to the calculated value 
of the manipulated steam flow so that part of the input signal was 
not correlated with the output product concentration, even though 
feedback control was used. By treating the dither signal as an 
input signal, the adequacy of the identified transfer function models 
for the closed loop system was determined by lack of correlation 
between the input and the final residuals, a. Dither signals with 
different characteristics were employed and are presented in Table 
3.13. Typical runs are shown in Figs. 3.14 to 3.18.. As expected, 
the output standard deviations for disturbances and dither signal 
addition under closed loop are smaller than those for open loop 
eondittons, (¢f. Tables 3.4 and 3.13): 

With the dither signal as an input and product concentration 
as output, closed loop process and noise models were estimated by 
proceeding in the same manner as for the open loop runs. Simulations 
with and without additional feed flow perturbations were carried out. 
Runs without feed flow disturbances were performed using three 
different controller constants for the product concentration steam 


loop but using the same level controller constants. The following 


control matrices were used: 
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TABLE 3.13 


CONDITIONS FOR CLOSED LOOP SIMULATION RUNS 


(Steam and feed flow disturbances, product 


concentration and level control) 


(1b /min ) 


SS-13 


SS-14 
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SS-16 


SSF-11 
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3.7.1 Model Estimation : Steam Disturbances 

Table 3.14 presents the models and parameter estimates for 
runs SS-11 to SS-16. For the runs SS-11 to SS-14, using 
controller matrix Kio 


closed loop transfer function models obtained are in very good 


with steam flow disturbances only, the 


agreement. The variations in parameter values for most of the 
parameters are less than 1 percent. However, attempts to carry out 
simultaneous transfer function and noise model estimation were 
unsuccessful for all these runs. Simultaneous estimation resulted 
in large chi-square statistics even after a large number of iterations 
(>40). Hence first the transfer function model and then the noise 
model for the residual series were calculated. The sum of squares 
of residuals after fitting the transfer function model alone was very 
small. 

For runs $s-11 and SS-12 the noise models are almost 
identical. It was not possible to fit reasonably any noise model 
with up to 5 parameters to the residuals from run SS-13. Further- 
more, the chi-square test on the cross correlations for this run 
indicates little confidence in the transfer function model even though 


the fit seemed reasonably good from the small residual sum of squares. 
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This was also the case for run SS-15. On the other hand, for 
run SS-16 with control matrix Ka reasonable process and 


noise models were obtained by simultaneous estimation. 


3.7.2 Comparison with Open Loop Models 


Open loop process models calculated from these closed loop 
models using eqn. (2.23) are presented in Table 3.15. Comparing 
these results with open loop process models calculated from open loop 
simulation data in Table 3.2, indicate that significantly different 
process model. orders are obtained from eT a Ci on under closed 
loop and open loop identifications. Autoregressive parameter S is 
nearly equal to the estimate obtained using open loop runs SS-4 and 


Lal 


SS-6. The estimate of 4 from the closed loop runs is almost equal 


2 
to the value obtained from run SS-6. Similarly, the estimates for tw 5) 
are comparable to those obtained from open loop runs. 

Although the closed loop runs have produced consistent 
parameter estimates (runs SS-11 to SS-15), the chi square test 


indicates very little confidence in the model orders. Hence the open 


loop models seem to be more reliable in this identification study. 


3.7.3 Model Estimation: Steam and Feed Flow Disturbances 

The parameter estimates for simulation runs with steam and 
feed flow disturbances are presented in Table 3.16. With the excep- 
tion of run SSF-11, the estimated process model parameters remain 
fairly constant with variations of less that 1 percent. The noise 


models in the five cases are second order with one degree of 
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FIGURE 3.17 CLOSED LOOP SIMULATION DATA FOR RUN SSF-14 
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FIGURE 3.18 CLOSED LOOP SIMULATION DATA FOR RUN SSF-15 
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difterencing. The noise model for. runs SS-11, 13, 15 and 17 
is (2, 1,1) and the parameters are within 10 percent of each 
other. Similarly for runs SSF-14 and 16 with noise model 


a “a 


(1, 1, 1), parameters b> 8) are nearly equal. 

For all the runs, the chi square statistics are fairly 
low indicating good fits were obtained. 

Table 3.17 presents the open loop process models calculated 
from the identified closed loop models for runs SSF-11 to SSF-17. 
As could be expected from the closed loop models, the parameter 
estimates here are almost equal, when Beet eed with the parameter 
estimates for open loop runs SSF-1 to SSF-4. In Table 3.6 the 
parameter estimates from the closed loop runs are within 5 ~- 10 
percent of the open loop values. 

Here, both identification techniques produce reliable models 


with adequately small values of Le and Q, indicating good overall 


tits, 


3.8 Closed Loop Identification Experiments 


Five runs, SF-11 to SF-15, were made using dither signals, 
feed flow disturbance sequences a feedback controller matrices. 
The. run conditions are given in Table 3.18 and typical plots are shown 
vrgih dis Medch) Newel RCE Gwe fee Be 


The various control matrices that were used are: 
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For run SF-15 an integral control matrix K, was used in addition 


= 
Co proportional control matrix, Ris 
0 OO 083562 
Ky = 07 20 0 
Oe) ae tO CeO da 0 
The proportional-integral control law used is 
i AT 
u(t) = Ky. x(t) + Kyte x (n)) > (3.11) 


As was observed in the closed loop simulations, the standard 
deviation of the product concentration was reduced when feedback 


control was used. 


3.8.1 Closed Loop Process Models 


Parameter estimates for the runs SF-11' to SF-15 are 
presented in Table 3.19. The process model with orders (r, s, b) 
a2 el et eaprovided good fits for runs’) SF-1) to SF=15° in spite 


of the different controllers used. For runs SF-11 to SF-13 with 
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control matrix K the o values were almost equal, being within 


1 percent of each other. The b5 values showed slightly larger 


“A n 


variations of approximately 5 percent. Parameters 05 and Wy 
obtained from runs SF-11 and SF-12 have comparable magnitudes but 
for runs SF-13 they differ by factors of 3-10. 

Table 3.19 indicates that although the fits are good for 
all runs, both with respect to the residual sum of squares and chi 
square tests, the model parameters are sensitive to change in process 


conditions and the control matrix. 


The closed loop study shows that models with transfer 


function model orders (2, 1, 1) have provided good fits for all runs 


but simulation runs SS-1l to SS-16. As discussed in honda. the 
164) are 10-15 percent higher in the case of models derived from 
experimental data, when the identical controller matrices are used. 


The w,} are higher for models estimated from simulation data. 


j 
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For closed loop study the simulation data with a steam dither signal only, 


produced models with different model orders and in many cases simultaneous 


fitting of transfer function and noise models was not possible. With 
this exception, other fits are generally good in both simulation 


and experimental studies. 


3.8.2 Comparison with Open Loop Models 


By using eqn. (2.23) open loop process models were 


calculated from the identified closed loop process models. These are 


presented in Table 3.20. The level control systems for runs SF-1 to 4 


(eigetaple oll) and? *SF=115912,"193Mi5s @ere*the same.e Comparison 
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a 


of these runs indicate that 8, from closed loop data is fairly 
close to the values from open loop data. Parameter 5, in runs 


SF-11 and SF-13 compare well with that for runs SF-1, 3, 4. 


W, and Wy values are of comparable magnitudes, but exhibit some 


variations, as was the case for the open loop runs. Ws is smaller 


“a 


than that obtained using open loop data, by about a factor of two. 

The small values of os in Tables 3.11 and 3.19 for these 
open loop and closed loop runs, indicate that in general good transfer 
function models were obtained. 

Since larger Q, values were obtained for the closed loop 
runs than for the open loop runs, this suggests that the noise models 


identified under closed loop conditions are less reliable. 


3.9 Conclusions 

Evaporator models calculated from both open loop and closed 
loop data demonstrate that the time series modelling technique 
presented by Box and Jenkins [1] can be successfully used to model 
real processes. The extremely small sum of squares of residuals 
indicated a very good fit to both experimental and simulation data. 

Process transfer function models obtained from simulation 
data containing "unknown" disturbances are fairly invariant whereas 
the noise models are sensitive to process conditions (e.g. type and 
magnitude of disturbances). For process models developed from 
experimental data, the estimated autoregressive parameters were 
constant (for given set of level control constants) but there were 
large variations in the moving average parameters. The fit of data 


is again very good, even in the presence of significant process 
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disturbances. 

The utility of the time series modelling technique was 
tested for processes under closed loop operating conditions. The 
resulting open loop process models calculated from the closed loop 
data were nearly the same as the open loop process models obtained 
from open loop data. Hence it can be concluded that time series 
analysis technique Broce to be a useful modelling technique for 
processes operating under open or closed loop conditions.. 

An overall comparison of the models, calculated from 
experimental and simulation data, indicates that for both open loop 
and closed loop conditions, comparable models were obtained. Higher 
6,4) values were obtained from the experimental study, whereas 
significantly higher tu,} estimates were calculated for the 
simulation data. No significant difference in noise model orders 
was observed. Similarly, with the exception of simulations using 
steam disturbances only, comparable and adequate fits were obtained 
using simultaneous transfer function and noise model Ettione 

The choice of input signal plays an important role in 
identification processes. Where the variance of input is constrained, 
use of first order autoregressive input signal with positive 
parameter 9 produced better excitation of the process and also 


reduces the width of confidence intervals on the parameter estimates. 
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CHAPTER FOUR 


EVALUATION OF MINIMUM MEAN SQUARE ERROR CONTROLLERS 


4,1 Introduction 

In this chapter controllers are designed using the evaporator 
models that were reported in Chapter Three. The noise models are used 
to forecast future disturbances and the transfer function models are 
used to take appropriate compensatory action. The resulting controllers 
are referred to as minimum mean square error (MMSE) controllers. 

In this chapter the performance of these MMSE controllers is 
evaluated and compared to the results obtained using conventional PI 
control or no control in the same disturbance environment. This 
comparison is carried out using simulated responses from the state space 
evaporator model and then using experimental data from the actual pilot 
plant evaporator. The robustness of the MMSE controllers is tested 


using deterministic disturbances such as_ step changes. 


4.2 Evaluation of MMSE Controllers : Simulation Study 


The MMSE controllers designed based on open loop identifica- 
tion runs SSF-1 to SSF-4 and closed loop runs SSF-11, 12 and 14, 
are presented in Table 4.1. The performance of these controllers was 
studied for stochastic disturbances and deterministic disturbances 
such as step changes in feed flow rate. 

If the noise model has a pole on the unit circle, then the 
resulting controller provides integral action [5] and thus eliminates 


offsets due to step disturbances. For this reason noise models with 
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TABLE 4.1 


MMSE CONTROLLERS DESIGNED USING MODELS FROM 


SIMULATION RUNS 


ieee a ee RN 


No. Run No. Controller, C(B) 


123.96 - 251.04 B + 160.28 B°~ 32.89 B° 


Pmt R6s4 B 0/4482 mor 0. 613608~ 


] SSF-1 


et ae 


108.68 - 233.86 B+ 161.29 BY - 35.91 B 
1 = 1.2630 B - 0.3785 BY + 0.6414 B> 
2 3 4 
DOD EABOMORBETEA (Oi) 2eBGl ml oad IBBre 5004 B 
1 - 1.8757 B + 0.2865 B¢ + 1.1196 B> - 0.5304 B° 
2 3 4 
218.68 - 492.33 B + 423.21 Be - 167.97 Bo + 26.53 B 
Hee 1084908 FOL 047 IBA. 12170 Be ONses1l Rs 
2 3 4 
144.23 B - 322.67 B + 256.53 B2 - 90.02 Bo + 11.95 B 
1 - 1.4301 B - 0.1748 B? + 0.9304 B> - 0.3255 B¢ 
2 3 i 
149.96 ~ 315.00 B + 225.22 B¢ - 68.35 BY + 7.76 B 


1 = 1.2792 B - 0.3014 B* + 0.7520 B> - 0.1714 BA 


113.27 — 224.80 B + 138.23 B* - 27.03 B> 
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FIGURE 4.1 SIMULATION RUN SFD-4 
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one degree of differencing were selected for controller design 
whenever such a noise model indicated good fit. Of the seven closed 
loop runs SSF-11 to SSF-17, three runs (SSF-11, 12 and 14) with 
different noise model orders were arbitrarily chosen for controller 
design. 

Simulation runs were made for four different first order 
autoregressive disturbances, SFD-1 to 4, and three different step 
changes in feed flow rate, SFS-1 to 3. The step changes consisted 
of a step increase at t = 5 anda step decrease of equal magnitude 
at t = 205. In all the runs the levels were controlled using 


proportional control matrix K The product concentration was 


1° 
under MMSE control with steam flow as the manipulated variable. The 
steam flow adjustments were constrained as -1.0 < u, < 0.5 where 
Lipa ies S)/S and $ is the normal steady state value. Two PI 
controller settings were used for the runs with PI controllers. The 
first set of controller constants, KP = -9.55 and. KI = -0.358, were 
obtained by King and McNeill [23] by experimentally tuning the product 
concentration control loop. A second set, KP = -40.6 and KI = -0.36 
was obtained from simulation data by tuning the PI controller to 
minimize the integral squared error (ISE) performance criterion for 


stochastic disturbances (oF =O. / Dy bp = 0.75). Runs with no 


concentration control are also presented for purposes of comparison. 


4.2.1 Performance Under Stochastic Disturbances 


The run conditions and the corresponding standard deviations 


obtained using 320 data points are presented in Table 4.2. Figures 
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FIGURE 4.2 SIMULATION RUN SFD-3 
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FIGURE 4.3 SIMULATION RUN SFD-4 


Th 


(38) 


> 
| _ 
+ re ; 
' i 
i 


‘ | ane) 
ih Wet , 


wh 


‘ : bide ‘ 
(CA Gad. 4 HA 


J 
ian? 
at | 


ASA NT Wa Vive dd 
Oba? (Ay? One 


oe ne So ee . SR Sete Ace Cee ean eens Deena 
ot N ae 
H W 
ra fe 
oar, us 
3 : pe Wag (22 
i 
¢ <a et 
a, sy 


Soe 
Se 
a 


NIWET NE Alvord 3 NIWG? NI MO WALLS 


G 
3 


ve 


240 


pha 2) 


TIMe IN MINS 


23 


“WIAD Jd Nf NOG PNP ND Ag 


PI (sim) (CONTROLLER 


FIGURE 4.4 SIMULATION RUN SFD-4 
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4.1 to 4.4 compare responses for open loop and PI controls and 
Fig. 4.5 and 4.6 show responses for the MMSE controller #2 for 
run conditions SFD-2 and SFD-4. Figures 4.7 and 4.8 present 
responses using MMSE controller #5, for the same set of random 
number sequences. 

It can be seen from Table 4.2 that as expected, both PI 
and MMSE controllers provide significant improvement over the open 
loop responses. Controller #4 based on run SSF-4 was found to be 
unstable under all run conditions, when no constraints were present 
on the manipulated variable. For this controller the response 
consisted of oscillations with slowly increasing amplitudes of the 
manipulated and controlled variable, thus indicating instability. When 
contraints were imposed on the manipulated variable the response 
in Fig.::4.9 was obtained. It can be clearly seen that the controller 
performance is not satisfactory since the manipulated variable 
frequently hits the constraints and the C2 response contains undesirable 
oscillations. The other three controllers based on open loop 
identifications performed well and provided better control than the 
two PI controllers. 

For the case of the MMSE controllers #5 to 7, based on 
closed loop identification runs, all the controllers performed ue 
well as the other MMSE controllers when stochastic disturbances 
were considered. As can be seen from Table 4.2, MMSE controllers 
produced standard deviations of errors, Sag? which are approximately 
one third of those obtained using PI (sim). The MMSE controllers 
based on closed loop identification runs produced more variations in 


steam flow as can be seen from Figs. 4.7 and 4.8. It can also be seen 
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from Table 4.2 that controller PI (sim) produced a standard 
deviation which is 60% lower than that for the PI (exp) controller 


tuned by King and McNeill [23]. 


4.2.2 Performance under Deterministic Disturbances 

The MMSE controllers were designed using models identified 
under stochastic disturbances similar to those used for runs SFD-1 
to SFD-4. Hence their improvement over PI control in the presence 
of stochastic disturbances is not surprising. To test the robustness 
of these controllers for other types of disturbances, runs with step 
changes in feed flow rate were made. Figures 4.10 to 4.14 show 
typical responses for condition SFS-3. From Table 4.2 it can be 
seen that the MMSE controllers performed better than PI(sim) and 
PI (exp) controllers for step changes as well. Here the reduction in the 
error standard deviation is about a factor of 3 for controllers based 
on runs SSF-1 and SSF-2 and as much as 20 times for the controller 
based on run SSF-3. The results clearly indicate that ability of 
the MMSE controllers to handle deterministic disturbances as well as 
stochastic disturbances. One of the MMSE controllers (#7) based on 
closed loop runs, produced higher standard deviation than that for 
PI controllers. Its response shows a steady increase in the 
deviation of product concentration from setpoint value, indicating 


unstable behaviour. 
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4.2.3 Performance under Setpoint Changes 

The previous study of MMSE controllers indicates that they 
perform better than PI controllers (in the sense of minimizing sum 
of squares of errors), when the problem is one of regulatory control. 
In order to further evaluate their performance, setpoint changes were 
also considered. In these runs a 1 percent setpoint change was intro- 
duced at the 30°F sampling interval and the standard deviations of 
errors produced by typical MMSE controller (#2) and PI (sim) controller 
were compared. Two types of responses were simulated for comparison: 
For the first run, feed disturbances were absent and in the second, a 


feed disturbance (o, = 0.75 1b/min, = 0.75) was present. Typical 


op 
results are presented in Figs. 4.15 and 4.16. 
The standard deviations of errors are shown below: 


Feed Concentration 
(1b/min) (% glycol) 


fo) 
FE C2 


PI (sim) MMSE #2 
0.08896 0.1002 
0.08813 0.09926 


Figure 4.16 ee He that MMSE controller #2 produced a 
stable response and from the above [9 values, it is obvious that 
the error standard deviation is only 10 percent larger than for well- 
tuned PI (sim) controllers. This study suggests that MMSE controllers 


can also be used successfully for setpoint changes. 
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4.3 Evaluation of MMSE Controllers : Experimental Study 

Minimum Mean Square Error controllers were designed, based 
on open loop experiments SI-L, 3 and 4. Since the noise models for 
these runs had a factor of (1 - B) in the denominator, these controllers 
include integral action. These MMSE controllers produced large input 
adjustments, hence some tuning was required before these controllers 
could be successfully implemented. In this section the required 
modifications and the controller performance are discussed and 
comparisons are made with the well tuned PI controller reported by 
King and McNeill [23]. As in the simulation study, both stochastic 


and deterministic load disturbances were employed. 


4.3.1 Implementation of MMSE Controllers 

The output signal of a MMSE controller is calculated as 
a weighted sum of present and past values of the error and the 
Manipulated variable. These controllers were implemented using the 
newly added 'Z-transform" controller option in the Direct Digital 
Control (DDC) package [24] available on the Department's IBM-1800 
computer [25]. This alogorithm provides storage for one gain term 
with ae weighting parameters for past values of errors and four for 
past values of the manipulated variable. Control was implemented 
by cascading the "Z-transform control" loop with steam flow control 
loop. Further details of the implementation are discussed in a 
research report by Kogekar and Berka [24]. 

The levels Wl and W2, were under proportional control 
using controller KI (eqn. 3.7). Level control and data acquisition 


for the experimental runs were carried out employing the multivariable 
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TABLE 4.3 
MMSE AND MODIFIED STOCHASTIC CONTROLLERS DESIGNED 


USING MODELS FROM EXPERIMENTAL OPEN-LOOP IDENTIFICATION 


Se 


MMSE Controller 
C(B) = K C (B) 
Modified Stochastic Controller, 


C'(B) = K' Cc, (B) 


jo. | mm tos) x 
epee 
pape 
1 = 1.9524 B + 1.1341 B* - 0.1779 B- 
% Tt jae 230.38 | 40.66 | et ee SS 
1 = 1.3874 B + 0.2932 B* + 0.0942 B 
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control and identification programs developed previously for the 


evaporator system [26]. 


4.3.2 The Modified Stochastic Controllers 

Table 4.3 presents the minimum mean square error controllers 
designed for open loop identification runs SF-l, 3 and 4. In each 
experimental run, the MMSE controller was initiated after the evaporator 
had attained the nominal steady state. However, even before the feed 
disturbances were introduced it was observed that for very small 
deviations in the product concentration (of the order of 0.01 - 0.02%) 
the steam flow rate manipulations were very large, reaching the upper 
and lower physical constraints quite frequently. Instead of reducing 
the output deviations from the setpoint value, the controllers produced 
an unstable condition. Hence it was necessary to modify the MMSE 
controllers in order to reduce the variations in the manipulated 
variable and avoid excessive corrective action for small deviations 
of the concentration. This was done by lowering the value of gain, 
K, by fifty percent or more, of the value obtained for MMSE controller. 
(The "gain" is defined by the expression in Table 4.3). The reduced 
"cain" values, K', are also presented in Table 4.3. These values 
were obtained by reducing the "gain'' in successive steps until 
adequate control, in absence of feed disturbances, was achieved. Thus 
for the modified stochastic controllers in Table 4.3, all the parameter 
values other than the "gain" K, have the same numerical values as 
the original MMSE controller parameters. 


The identified evaporator models have a root near the unit 


circle (e.g. for run SF-4 the roots are 1.0185. 1.393). This means 
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a change in the manipulated variable has little immediate effect on 
the controlled variable [1, p. 472]. Therefore the MMSE controller 
calls for excessive input manipulations. In conditions such as this, 
modifying the controller is recommended [1, 11, 12]. The constants in 


the numerator polynomial w(B) have values of the order of lta to 


107" and these values change significantly with changes in process 
conditions. The 95 percent confidence limits on these parameters 

are very wide (c.f. Appendix D). The effect is that the numerator of 

the controller has parameters of the order of 10° and the "gain" 

K is subject to large variations with wide 95 percent confidence limits. 


Hence, in order to have a stable stochastic controller a reduction 


in the value of "gain" term was attempted. 


4.3.3 Discussion of the Results 

The performance of stochastic controllers obtained from 
various identification runs was evaluated for two stochastic CEeirst 
order autoregressive) disturbances and for an approximately 20 percent 
step decrease in feed flow rate from 5.5 1lb/min to 4.5 1b/min. The 
various run conditions and the standard deviations of feed flow, steam 
and product concentration are presented in Table 4.4 and Figs. 4.19 - 
4.30. 

During the experimental control runs, periodic fluctuations 
in the building vacuum supply were observed. In order to ensure that 
these fluctuations were external to the evaporator apparatus, the 
building vacuum supply was also recorded in the Department's Instrument 
Shop and is presented in Fig. 4.17. Approximately every 100 minutes 


there was a sudden drop in the vacuum from -21.3 inch Hg to ~-20).3 inch, 
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FIGURE 4.17 CHART RECORD OF BUILDING VACUUM SUPPLY PRESSURE. 
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18 EFFECT OF VACUUM DISTURBANCES ON PRODUCT CONCENTRATION 


(concentration) 
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then the vacuum slowly increased back to -21.3 inch during the 
next 100 minutes. In the evaporator system, the vacuum passes through 
two vacuum damper tanks in series before reaching the Separator. The 
separator pressure is controlled at a value of -18.75 inch Hg by 
manipulating a slip stream of air entering the separator. Even under 
these conditions the sudden drop in vacuum supply caused a pressure 
increase of approximately 0.25 - 0.50 inch Hg in the Separator. Since 
the product concentration is very sensitive to variations in second 
effect pressure, P2, these periodic pressure increases caused sharp 
decreases in the product concentration as well as decreases in the 
separator ieee The vacuum behaviour is shown in Fig. 4.17 and its 
effect on concentration is illustrated in Fig. 4.18 by presenting the 
analog plots of product concentration and P2. (The problem 
with the building supply pressure was not present during the identifi- 
cation runs described in Chapter Three). 

During the control runs these sudden decreases in concentra- 
tion caused corresponding sudden increases in steam flow rate. As a 
result, oscillations in steam and concentration resulted and the 
standard deviation of errors increased. It was. observed that the 
stochastic controllers performed poorly whenever the concentration 
suddenly deviated from setpoint value by more than 5 percent glycol. 
For anit changes, small offsets resulted whereas for larger 
deviations, unstable behaviour was observed. By contrast, the PI 
controller provided satisfactory control (c.f. Figs. 4.23 and 4.27). 

In summary, the control would have been much better if the 


vacuum had remained constant and the stochastic controllers did not 


In Figs. 4.19 to 4.30 arrows indicate drop tn vacuwn. 
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prove very useful for handling these disturbances which tended to 
resemble setpoint changes. 

Table 4.4 shows that the error standard deviation for PI 
or stochastic control is 4 - 7 times lower than the values for open 
loop (or no control). Compared to PI controllers, modified stochastic 
controllers produced a 10 ~ 50 percent reduction in standard deviation 
for the same flow disturbances. The mean value of the product concentra- 
tion is almost the same as the setpoint indicating no offset for 
stochastic disturbances. 

Although the standard deviation of. product concentration has 
decreased when stochastic control is used the steam standard deviation 
has increased by approximately 50 percent, with the exception of run 
FD-23. 

Runs FS-11 to FS-14 were carried out for a -20 percent 
step decrease in feed flow rate. The standard deviation was calculated 
starting form the point at which step change was introduced. Here 
also the standard deviations of the errors are smaller than those 
obtained using PI control by 10 - 50 percent, in case of runs FS-12 
and FS~-13. The steam standard deviations were nearly equal for runs 
FS-11 and FS-13 and were 25 percent more for run FS-12. Here the 
deviations of the mean from the setpoint are more than for stochastic 
disturbances. | 

It can be seen from Figs. 4.18 to 4.20 that the final 
steady state values are different from the setpoint or the initial 
steady state value. But the deviations did not generally begin from 


the time when the step change in feed flow is introduced as they were 


caused by the sudden changes in vacuum. 
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When the three controllers were designed, noise models 
with one degree of differencing were used to provide integral action 
in the controllers and remove offsets. But since the process model 
has a root very near the unit circle (1.01) the effect of the (1 - B) 
factor providing integral action is almost nullified and the result 


is a controller without integral action. 
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Identification 

The MMSE controllers designed using closed loop identification 
runs SF~l1 to SF-14 are presented in Table 4.5. In comparison 
with the controllers in Table 4.3, these controllers have higher K 
values. Various runs were attempted using these controllers and 
successively decreasing the value of K to ‘19 but for all these 
runs, the stochastic controllers failed to maintain the steady state 
concentration even before the feed disturbances were introduced (i.e. 
unstable responses were obtained). 

The roots of the closed loop system in Fig. 2.4 when MMSE 
Sehetor tees #4 - 6 were used, were studied. With the exception of MMSE 
#4, when each of the two other MMSE controllers was considered, the roots 
of the closed loop system were outside the unit circle indicating a stable 
closed loop system. When MMSE #4 was used, it had a root (0.7033) 
inside the unit circle indicating unstable behaviour. 

Simulations carried out using the MMSE controllers in Table 
4.5 produced oscillatory and unstable responses, both for constrained and 


unconstrained manipulated variables. The oscillations persisted even 


after reducing K by fifty percent or more. 
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TABLE 4.5 
STOCHASTIC CONTROLLERS DESIGNED FROM EXPERIMENTAL 


CLOSED LOOP IDENTIFICATION RUNS 


(Controller C(B) = K C (B)) 


No. | Run No. K C (B) 


— 


1 - 1.920 B + 1.2698 B- - 0.3541 B> + 0.0939 B4 
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4.4 Conclusions 

From the study of MMSE controller performance, it was 
concluded that in simulation study the MMSE controllers, in all but 
one case, handled the problem of regulatory control under stochastic 
and deterministic disturbances very well. Only MMSE controller #4 
produced unstable responses It was also demonstrated that performance 
comparable to PI controllers was obtained when setpoint changes were 
introduced. On the other hand, the controllers based on the 
experimental identification runs needed some tuning before they could 
be successfully used. In spite of vacuum disturbances they 
performed, in most cases, better than PI controllers. Severe vacuum 
disturbances led to different final values than the setpoint value. 
Because of small integral action, these seem like an offset for short 
runs. Observation of the responses indicates that the concentrations 
are slowly approaching the setpoint value. 

For systems with low transfer function gains and large ie 
constants, tuning or modifying the MMSE controller may be required to 
provide smaller input variations. 

Although the controllers based on closed loop identification 
performed well in the simulation study, they failed to 
produce stable responses in the experimental study. This indicated less 
reliable process model identification under closed loop conditions from 


the point of view of controller performance. 
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CHAPTER FIVE 


CONCLUSIONS AND RECOMMENDATIONS 


The effectiveness of the time series modelling technique 
proposed by Box and Jenkins [1] has been investigated using simulated 
and experimental data for a pilot scale evaporator. The results 
clearly demonstrate that this technique can be successfully employed 
to obtain simple, dynamic stochastic models for a high order, inter- 
acting process such as the double effect evaporator. Identification 
using open loop and closed loop data, produced open loop process models 
that are in good agreement and exhibit similar trends. Comparable sum 
of squares of residuals and chi square statistics indicate that both 
identification methods are equally efficient. 

Transfer function models with model orders of (2, 1, 1) 
proved adequate to model the process under various types of disturbances. 
Autoregressive parameter estimates 65) were fairly insensitive to 
process conditions and also to the choice of noise model orders, whereas 
the moving average parameters, a) were sensitive to changes in 
process conditions and to the choice of noise models. Different noise 
models were required to model different disturbances. 

MMSE controllers designed using these models produced stable 
responses in the simulation study, except for one controller. However, 
the experimental implementation of MMSE controllers caused large 
adjustments in the input. Due to physical constraints on the input, 
tuning the controller to reduce the input variance was necessary to 
achieve stable responses. The controllers based on closed loop simula~- 


tion data performed well in simulations, but controllers designed 
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using models identified from closed loop experimental data were found 
to be unstable even after large gain reduction. 

All stochastic controllers exhibiting stable responses 
produced smaller sum square of errors for both stochastic and determin- 
istic disturbances than the well tuned PI controllers. However the 
PI controllers handled vacuum disturbances better. In the simulation 
study, MMSE control was successfully employed for setpoint changes. 
This indicates the robustness of these stochastic controllers. 

It can be concluded that the time series modelling technique 
can be employed to obtain simple models for processes. Robust 
stochastic controllers can be designed from these models, but in some 
cases the controllers produce highly oscillatory, unstable responses. 

By comparing the error standard deviations using PI control, 
(2 parameters) and MMSE control (5~7 parameters) in the experimental 
study, the small improvement in control does not seem to justify, 
from point of view of industrial applications, the use of a larger 


number (5 - 7) of parameters. 


5.1 Recommendations for Future Work 

1. For the double effect evaporator application, by treating 
feed flow as a measured disturbance, feedforward - feedback control 
can be implemented. 

2. The single-input single-output analysis used in this 


study can be extended to multivariable systems, taking into account the 


process interactions. 
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3. As some of the MMSE controllers produced unstable or 
highly oscillatory response, a technique to determine the stability 
of the closed loop system may be developed. 

4, A method of reducing the number of controller constants, 
without increasing the error variance significantly, needs to be 
developed to achieve wider applicability. 

5. It was found for the application of MMSE controllers to 
the evaporator that the input variations required were large. In 
order to reduce the input adjustments a controller that minimises the 
quadratic cost function E fe + AXE} may be designed using the 


evaporator models. 
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NOMENCLATURE 
(a) Alphabetic 
a White noise 
B Backward shift operator 
b Dead time as whole multiples of sampling interval 
Bl First effect bottoms flow rate 
B2 Second effect bottoms flow rate 
Cc Output coefficient matrix 
Cl First effect concentration 
C2 Second effect concentration 
CF Feed concentration 
D Dither signal 
d Degree of differencing 
d Disturbance vector 
E Expectation operator 
e Error 
F Feed flow rate 
Hl First effect enthalpy 
HF Feed enthalpy 
K Gain 
K Gain matrix 
Ny Process noise 
01 Overhead first effect 
02 Overhead second effect 
Pp Order of autoregressive part of noise model 


Pl First effect pressure 
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NOMENCLATURE (continued) 


PZ Second effect pressure 

Q Chi-square statistic for autocorrelation function 
q Order of moving average part of the noise model 

r Order of denominator polynomial of the transfer 


function model 


aay Cross correlation function between series 
bg and WY, } 

8 Steam flow rate 

S Chi-square statistic for cross correlation function 

ath Temperature of first effect 

TZ Temperature of second effect 

TF Feed temperature 

is Sampling instant 

u Control vector 

V Impulse weight polynomial 

Wl First effect holdup 

W2 Second effect holdup 

xX Input variable 

24 State vector 

nM Output variable 

a Output vector 

(b) Greek 

a Prewhitened input 


B Transformed output 
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NOMENCLATURE (continued) 


X Chi-square statistic 

6 Polynomial 

) Autoregressive polynomial in noise model 
i) Partial autocorrelation 

Y Autocorrelation function 

A Polynomial 

mT Mean value 

n Polynomial 

0 Correlation function 

x Summation operator 

Oo Standard deviation 

Wy Polynomial 

8 Moving average polynomial in noise model 


Disturbance coefficient matrix 


j@ 


(c) Subscripts 


a ae element 
th 
G t sampling instant 
x Pertaining to variable x, (i.e. oe standard deviation 
of x) 
ss Steady state 
= Vector 
Matrix 


aa Between series to} and la, } 
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NOMENCLATURE (continued) 


(d) Superscripts 


=] 


A 


Matrix inverse 
Estimated value 
Normalized value 


Matrix transpose 


(e) Abbreviations 


iefelp es, 
Cree 


CC 


C.I. 


MMSE 


Dea. Geils 


Auto correlation function 
Cross correlation function 
Concentration controller 
Concentration recorder 
Confidence interval 

Flow controller 

Flow recorder 

Integral constant 
Proportional constant 
Level controller 

Level recorder 

Minimum mean square error 


Proportional plus integral 


Partial auto correlation function 
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APPENDIX A 


THE EVAPORATOR MODEL 


The fifth-order discrete state-space evaporator model 


calculated by Wilson [20] is represented as: 
(to) =o Ct) a Ce) 8d (b) (Aj) 
HAGE STE Res, (A. 2) 


The elements of the vectors x, u, d, y are defined 


as normalized perturbation variables, e.g. 


Wl - Wi 
Souci, 


Wl 
ss 


a kee! 


Where Wie is the normal steady state value of Wl. The vectors 


By) U; ad and y are defined as follows: 


State Vector, x Normal Steady State Value 
x = (Wl, Hey whale We) C2 

Wl First effect holdup - 45.5.1b 

Cl First effect concentration eI LE LVCOL 

Hl First effect enthalpy 189.2 BTU/1b 

W2 Second effect holdup 41.3 1b 

C2 Second effect concentration 1OCLLZ eiycol 
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Control Vector, u 

fieeoehsw el. 1B? | 

S Steam flow 

Bl First effect bottoms flow 


B2 Second effect bottoms flow 


Disturbance Vector, d 


d'=(F, cr, HF] 


F Feed Flow 
CF Feed concentration 


HF Feed enthalpy 


Qutput Vector, y 


y = (Wl, w2, c2] 


The coefficient matrices of the discrete-time model, with 


2 OF ib min 
3.485 1b /min 


1.581 1b /min 
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S27, elLycoL 


156.9 BIU/1b 


a 64 second time base, are shown in Table A.1. 
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APPENDIX B 
CALCULATION OF OPEN LOOP PROCESS MODEL FROM A 


MODEL IDENTIFIED UNDER CLOSED LOOP CONDITIONS 


Figure 2.3 and eqn. (2.23) indicate that the closed loop model relat- 


ing dither signal DY and output, Gs is the identified model: 


CpG (B) Deer E(B irae (B. 1) 
where 
eT Virani ek) See 
BO Srey are) Abe 
2 v(B) 
Ee) Sahl a—tC(BYAV) | ee? 
and 


C(B) = controller transfer function. 

V(B) = open loop process transfer function model. 
v(B) = open loop noise model. 

G(B) = closed loop process transfer function model. 
H(B) = noise model obtained under closed loop 


conditions. 


Equation (B.2) can be rearranged as, 


atic 
V(B)= TF C(B) GB) ee) 
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For closed loop experiment run # SF~13 in Table 3.19 


~3 
(1.08 + 0. 

G(py = (1:08 + 0.378 B) x 10 B (rs) 
1 - 1.697 B + 0.755 B 


The proportional controller gain used during closed loop operation 


was K = -2.5. Then substituting eqns. (B.5) and (B.6) into (B.4) 


gives: 
a3 
V(B) = CleU8pte0.57505 ub ex = (B.6) 
1 - 1.70 B + 0.754 B 
and 


Ie-nVCB) CCB) = 


Wewie700 Be 60.754 8-4 2.5 (1.08 Be 01378 B-) on 10. 


5 (Bow) 
bee 032 b 02/5458 


Since the numerator and denominator of (1 - V(B) C(B)) 


are almost identical, 
tee Vin) c(Ry =" 1 (B.8) 


hence 
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v(B) = H(B) = (79,4821 B - 0.2106 B) (1 - B) 


(B.9) 


The complete open loop process model is: 
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APPENDIX C 


MMSE CONTROLLER DESIGN 


The process model identified for simulation run SSF-1 
is presented in Table 3.4. Using the procedure outlined in 82.3 a 
minimum mean square error feedback controller is designed as follows. 
The process model has orders (2, 1, 1) for the transfer 
function part and orders (1, 1, 1) for the noise model. 


._ 0.01363 + 0.009609 B 


Y¥ x +N Ge 1) 


SP al ALE te Ma ae). Te 
and 
ee ACL 0.1808 B) 

Free mee0.870408 (iB yieat Ce 
Therefore 

Nee “2 bert ERS ee 1.6896 - 0.8704 B = a, (c.3) 

(1 = 1.87048 + 0.8704 B ) 
hence 
: - 0.8704 B 
C1=2 1, 87 045Re 0.38/04 ena) 

and 


A\(B) = 1 Geley. 
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Using eqn. (2.30) the MMSE controller obtained is then, 


75 (B)_n(B) 
NS) a (BYE) 


Peierls Be 0. S151 EB Ay nC1 e686 <1 0.87 00NB) 
(0.01363 + 0.009609 B) (1 - 1.8704 + 0.8704 B*) 


e303 6o0N 02s be 12930 8] 0. 2654) 
(l= 1.1654 B ~ 0.4482 Bo + 0.6136 B) 


(C.6) 
The control law for this controller is, 


=F in a - ae : 
K.4 = 7123.96 (e, - 2.025 e, | + 1.293 e,_, ~ 0.2654 e_4) 


an (1.1654 xX, + 0.4482 x5 = 0.6136 ia? 
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APPENDIX D 


RESULTS FROM PARAMETER ESTIMATION 
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(eisOnie) “Wla77? =Ove777) -0,03175 O.Olvae 1.0110 == 0.3525 
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St et Sam met FRY SeMd em mes ret eh WON HF0H Fah Se yee SOT St me trey Wnt bm pine Sfmt eT ge fry Sed Mm NY Mm ret Mh fee SPT WN BY Yom ts pau nem Pant Pent meee bee mest mene mt HE vs Sh tum te Pain mm ntry et a oa phy Rew Fame OR} Sem bem FaLH re SO rr men pe te ris 4 mg Km Meee Te rut at a nm i eS ee 
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CamOm tan) 1-628 0.6562 0.02487 0.003887 O.2748 0.28991 -0,3489 
1,628 —-0.6392 0600474 -O0.01707 0.45607 0.0605 0.5997 


‘wm seen sare vo Frm dAdo POR Y00 Pony re OP sV0m fem simy ROmR FOmD rest p/m Eman FvEN c¥ve Meat WOnD MET J+15 vamp fb) Sour rose FYM rte meat bnry FLAS Tur BPVE ntey rtnN Aber foun stm min Pee fee mm Rd Ir Heit Meat veh M4 mine Ofen arne eate wary tems S9m8 Me creme 8A Ange ome Mis 4A ym Some Reh ner Sg me oem RG H4 Hom nes tw 


1 Process transfer function model orders, (r,s, b ). 


2 Noise model orders, ( p, d, q ) 


** 95 percent confidence limits for the parameter estimates. 
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aoe na 6, 
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